Temperature dependence of critical velocities and long time decay of supercurrents 
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We show that a microscopic theory of superfluidity, based on the properties of the many-body 
spectrum of a superfluid, can explain naturally the temperature dependence of critical velocities 
and the long time decay of supercurrents. 

PACS numbers: 67.25.-k 



Landau gave the first understanding of persistent cur- 
rents in 4 He by relating this phenomenon with properties 
of the quasiparticle spectrum of the system [l| . Bogoli- 
ubov explained the linear dispersion of the quasiparticle 
spectrum at small momenta with the assumption of BEC 
fraction in a superfluid 0] , and thus supported Landau's 
theory. This microscopic picture of superfluid is, how- 
ever, not invoked in understanding some later important 
observations of superfluids, such as, *A) the temperature 
dependence of critical velocities 0,11], and *B) the long 
time decay of supercurrents 0] . It is natural to wish that 
a theory of superfluidity can account for these observa- 
tions besides explaining superfluidity, for the purpose of 
a unified picture to describe superfluidity and its related 
properties. In this paper, we show that a rather simple 
microscopic theory of superfluidity can accomplish such 
a wish. We also make comparisons between this theory 
and lordanskii-Langer-Fisher (ILF) theory 0], which is 
formulated to understand the temperature dependence of 
critical velocities and in which multiple assumptions of 
vortex rings are essential. 



This microscopic theory relates superfluidity with the 
properties of the many-body spectrum of a superfluid. 
The main conclusions of the theory are the following: 
i) The many-body dispersion spectrum of a superfluid 
E = E(P), where E is the lowest many-body eigen- 
energy at given momentum P, is not a monotonic func- 
tion of P; there exist energy barriers in the many-body 
dispersion spectrum which separate and prevent some 
current-carrying states from decaying 0, H, 0] (see Fig. 
[T]). ii) the existence of the energy barriers is due to Bose 
exchange symmetry 0]. Hi) The height of barriers will 
decrease with increasing momentum (velocity) of the cor- 
responding supercurrents; and beyond a certain velocity, 
the barriers disappear and the system dissipates its mo- 
mentum like a normal system p], 0] ■ We shall show that 
with this knowledge of the many-body spectrum, the ob- 
servations (*A, -kB) can be naturally explained. We first 
use the spectrum of a one-dimensional superfluid to il- 
lustrate these observations, and later discuss higher di- 
mensional superfluids of which the gross structures of the 
spectra are also the same as what specified above. 



We consider N Bose particles moving in a ring with a 




FIG. 1: A many-body dispersion spectrum of a ID Bose 
system. The barrier height decreases with increasing of the 
momentum at the corresponding minimum. The height of the 
barrier for the k-th minimum next to the ground state regime 
is Afc = Ao(l — k/(v c MR/h)) 2 with the momentum at the 
minimum is P = kNh/R. The minima disappear beyond a 
certain momentum. 



radial size of R. The Hamiltonian has the form of 



H 



N 

E 



d 2 



2MB. 2 86? 



+ 



N 

E 



(i) 



where 0j is the angular coordinate of the zth particle, 
M is the mass of a particle and V is the repulsive in- 
terparticle interaction. V can be either short-range or 
zero-range. The many-body spectrum can be classified 
by their angular momenta L (in following we replace L 
by a corresponding momentum P = L/R, considering 
the equivalent system moving in a straight line with pe- 
riodic boundary condition). Due to Galileo invariance, 
the full spectrum can be obtained if the part at momen- 
tum regime < P < Nh/R is known 0,0]. Specifically, 
if one denotes the n-th level at momentum P by e n (P), 
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then 

e n (P + kNh/R) = e n (P) + ((P + kNh/R) 2 -P 2 )/2NM, 

(2) 

where k is an integer. 

The many-body dispersion spectrum, as a function of 
P at regime < P < NK/ R, relative to the ground state 
energy, can be written in a form of 

E(P) = E ex (P) + HP/2MR, (3) 

where HP/2MR is a (trivial) part of kinetic energy, and 
where E ex (P) is symmetric (in E — P plane) to the 
line P — NH/2R and reaches its maximum at P = 
NH/2R. Numerical results @, [l(| suggest that E EX {P) 
is parabolic, i.e., E ex {P) = v c (Nh/R-P)P/(Nh/R), we 
shall use this form of E ex (P) for discussions. 

With the knowledge of E ex (P), one can obtain the full 
many-body dispersion spectrum and realize the follow- 
ing feature of this spectrum, i) for a positive integer 
k < v c MR/h, there is an energy barrier in the dispersion 
spectrum within the momentum regime (fc — l)Nh/R < 
P < kNh/R. ii) the height of this barrier is A& = 
A (l - k/(v c MR/h)) 2 , where A = v c Nh/4R (see Fig. 
P). 

Each valley of the dispersion spectrum can 'trap' 
metastable states of the system given that the temper- 
ature is lower enough, for a certain valley, we consider 
the quasi thermal equilibrium distribution of many-body 
levels in the region plotted in Fig. [5J then the possibility 
for the system staying at the levels lower than the energy 
barrier A is that 

e -(e s -e )/fc B T 

7= g e -( ei -o/* B r W 

Where i,j refers to levels in the region, ks is the Boltz- 
mann constant, and e is the lowest eigen-energy at the 
valley, i.e., the local minimum of the dispersion spectrum. 

At any small but finite T, 7 is smaller than unity. It is 
with some possibilities that some states, which are near 
the edge of the valley regime and close to the next lower 
valley (thick bars in Fig. [2]), can be reached, and sub- 
sequently these states can decay with ease to the states 
in the next lower valley, (i.e., they are kind of 'doorway' 
states). This process is then repeated and eventually 
the system transfers from this valley to the next valley. 
When 7 is the close to one, this transferring process is 
slow, which, nevertheless, leads to the decay of the su- 
percurrent at the large time scale 0, [3] . 

When 7 exceed a certain value j c close to unity, one 
can ignore the long time slow decay of supercurrents if 
only relatively short time scale is concerned, the value 
of 7 depends on the many-body spectra, T and A&. If 
one roughly estimates that 7 = "f c at T — cA^ at the 
k-th valley, where c is a constant, one then finds the 
highest possible velocity of supercurrents at T is given 
by v c (T) = v c (l — y/T/Tx) where T\ is the transition 
temperature [Xlj | . One shall note that the local spectra 
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FIG. 2: With a sufficient small temperature, a system can stay 
in a valley for a long time. The many-body levels (bars) in 
the valley region, which is bordered by the valley and roughly 
bordered by the thick dashed line, will reach a local quasi 
thermal equilibrium and the probability at a level with energy 
e„ is determined by the usual law, p oc e~ e "/ kBT . With a 
finite T, the states (thick bars) near the next lower valley are 
partially accessed and subsequently decay to the next valley, 
thus causing a leak of the system from this valley to the lower 
one. 

at different valleys are not the same, for example, the 
shapes of the valleys become more asymmetric with val- 
leys further away from the ground state regime. For this 
reason, the rough temperature-critical-velocity relation 
could be rendered to take a form of 

v c {T) = v c (l - ^ a (T)T/a(T x )T x ) (5) 

where a(T) is a slowly varying function of T. 

We thus illustrated the theoretical pictures of two ob- 
servations (*A, -kB) using the many-body spectrum. One 
could note that the pictures are direct and unavoidable. 
A spectrum of a system is a fundamental property of 
the system, and one can map physical processes to the 
transferring processes among the many-body levels of the 
system. 

Previously, ILF theory is constructed to explain the 
temperature dependence of critical velocities. It is in- 
teresting to compare our theory with ILF theory. One 
can see the following agreements between them: a) that 
there are energy barriers to prevent the decay of cur- 
rents at short time scale; b) the energy barriers can be 
'overcame' by thermal excitations, which leads to long 
time decay of supercurrents; c) the energy barrier is a 
decreasing function of the velocity of the supercurrcnt. 

The differences between two theories are the follow- 
ing: a) ILF involves multiple assumptions of vortex rings, 
such as their sizes, their dynamics, the creation and an- 
nihilation of vortex rings. Our theory suggests that these 
assumptions are not necessary [12j. b) ILF theory deter- 
mines the heights of the energy barriers, using energet- 
ics of vortices, and suggests the heights decrease linearly 
as the function of the_ velocity of the supercurrent, i.e., 
A(v) = A — p v 13]. these 'conclusions' lead to in- 
adequacies of the theory in its quantitative description 
of some systems [3]. They also lead to inadequacies in 
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FIG. 3: Possible scenarios of dispersion structure at the mo- 
mentum regime < P < Nh/R. E ax {P) = E{P)~ hP/2MR, 
where E(P) is the dispersion spectrum. (Upper panel) If cer- 
tain conditions on boundary and on N value are satisfied, 
there could be local minima at h/a, 2h/a,... (h = 2-kK). 
(Lower panel) In general cases, only the primary type of min- 
ima plays an important role in determining superfiuidic prop- 
erties. 

quantitative description at the temperature regime far 
below the transition point 4]. Within our theory, the 
height of energy barriers are naturally determined by the 
many-body spectra for which all low-lying eigenstates in- 
cluding the many-body dispersion states are relevant. 



We shall discuss higher dimensional cases. We consider 
for example a system of N particles in a tube with pe- 
riodic boundary condition. The section area a of a tube 
could be much smaller than R 2 , where R is the length of 
tube divided by 27r, but is order(s) of magnitude larger 
than a 2 , where a is the average interparticle distance. 
Again, with the knowledge of dispersion s pec tra in the 
momentum regime < P < Nh/R [7j, |9|, |l5j, one can 
derive the full dispersion spectra, particularly all local 
minima and the energy barrier height associated with 
each minimum, due to the Galileo invariance. 

Beside the primary type of minimum at P = Nh/R, 
with the corresponding many-body state in principle al- 
lows a fraction of BEC, in [l5| we also find some local 
minima at P = h/a, 2/i/a, the existence of these local 
minima has sensitive dependences on boundary condi- 
tions and on the value of N (see Fig. [3]). It also requires 
strong interaction. Once these conditions are not satis- 
fied, that type of minima at P = h/a, 2h/a, ... may not 
exist or be very shallow, and only the primary type of 
minimum plays an essential role in superfluidity. 

In conclusion, we illustrate that a microscopic theory 
of superfluidity, based on the properties of a superfluid's 
many-body spectrum, explains naturally the tempera- 
ture dependence of critical velocities and the long time 
decay of supercurrents. We emphasis that these pictures 
of a super-fluid are direct and unavoidable, due to that 
physical processes can be viewed as changes of occupation 
probability of each eigen-level and the transfers among 
them. 
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